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We investigate the electrical conductivity (a) of quark matter via the Kubo formula at finite tem- 
perature and zero quark density (T ^ 0, ji = 0) in the presence of an external strong magnetic field. 
For this purpose, we employ the dilute instanton-liquid model, taking into account its temperature 
modification with the trivial- holo no my caloron distribution. By doing that, the momentum and 
temperature dependences for the effective quark mass and model renormalization scale are carefully 
evaluated. From the numerical results, it turns out that a ~ (0.02 ~ 0.15) fm -1 for T — (0 ~ 400) 
MeV with the relaxation time r — (0.3 ~ 0.9) fm. In addition, we also parameterize the electrical 
conductivity as a/T « (0.46, 0.77, 1.08, 1.39) C EM for r = (0.3, 0.5, 0.7, 0.9) fm, respectively. These 
results are well compatible with other theoretical estimations, including those from the lattice QCD 
simulations. It also turns out that the external magnetic field plays only a minor role for a even 
for the very strong one Bo ~ nn\ x 10 and becomes relatively effective for T < 200 MeV. More- 
over, we compute the soft photon emission rate from the quark-gluon plasma, using the electrical 
conductivity calculated. 
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JL QCD at extreme conditions has been one of the most celebrated subjects in high-energy physics for decades. The 

^J" 1 observation of the strong magnetic field Bq ~ in the non-central heavy-ion collision at Relativistic Heavy Ion 
Collider (RHIC) of BNL [TH1] triggered abundant research works, such as the chiral magnetic effect (CME) [SHE], 
chiral magnetic wave (CMW) [9] [10], QCD phase structure in the presence of magnetic fields, related to the magnetic 
catalysis [TTHT3] , electromagnetic properties of QCD matter [T4lfT9] , and so on. In addition to the magnetic field, the 
transport coefficients for the hot and/or dense matter is of great importance in general as well, since they characterize 
the physical properties of the matter. In general, those coefficient can be studied by the Kubo formula [2QU22] . 
f^. Among them, in the present work, we want to investigate the electrical conductivity (a), which corresponds to the 
vector- current correlation (VCC) and its spectral functions [18], for quark matter at finite temperature and zero quark 
density (T ^ 0,/i = 0). This condition resembles the present heavy-ion collision experimental conditions, such as 

— RHIC and Large Hadron Collider (LHC) of CERN, in the presence of the external magnetic field. We note that the 
conductivity effects much on the magnetic field in terms of the relaxation time [14] [19] . The spectral functions are 
deeply related to the thermal dilepton production, and the production cross section contains much information on 
the hadron properties, produced in quark-gluon plasma (QGP) [23] . 

^ The electrical conductivity was investigated in the hot phase of the QCD plasma and extracted from a quenched 

SU(3) lattice QCD (LQCD) simulation of the Euclidean time VCC for 1.5 < T/T c < 3, resulting in a/T « 7C EM , 
where Cem = XI/ e / an d e f denotes the electric charge for the quark with a flavor /, in Ref. [22]. Similarly, the 
electrical conductivity was studied with quenched SU(3) LQCD with the maximum entropy method (MEM) for VCC 
in the deconfined phase T ~ 1.45 T c in Ref. [18 , in which the dilepton production rate was studied as well. In that 
work, the electrical conductivity turned out to be 1/3 < cf/(CemT) < 1. Note that the estimated value for a is about 
one order larger in Ref. [22] in comparison to those from Ref. [18 . In Ref. [19], VCC in quenched SU(2) LQCD 
simulation was investigated with a chirally invariant lattice Dirac operator with a constant external magnetic field. 
It was observed that a± and cry, which are perpendicular and parallel to the external magnetic field, respectively, 
are almost the same and insensitive to the external magnetic field for T > T c , whereas there appear considerable 
differences between them at T ~ with respect to the magnetic field. The electrical conductivity was also estimated 
as cr ± w a\\ « cr = (0.076 ± 0.010) fm -1 for T « 350 MeV with T = 1.12 T c . Beside the LQCD simulations, using 
the Green-function method, the electrical conductivity was explored at finite T and \i\ T « 100 MeV and \i ~ 400 
MeV [14 . At specific choices for the model parameters, the electrical conductivity was resulted in a w 0.04 fm -1 , 



* E-mail: sinam@kias.re.kr 



Typeset by REVTeX 



2 



which is similar to those from Refs. p~8j[T9]. Interestingly, they found that the electrical conductivity is very insensitive 
to the external magnetic field. 

Considering the rather unsettled situation for the electrical conductivity for hot and dense quark matter as men- 
tioned above, in the present work, we employ the dilute instanton-liquid model to study it [24] [25]. The model is 
characterized by the two phenomenological (instanton) parameters: average inter- (anti)instanton distance R ~ 1 fm 
and (anti) instanton size p ~ 1/3 fm. In the instanton ensemble, representing the nonperturbative QCD vacuum, 
the quarks are delocalized by flipping their helicities. As a result, the quarks acquire their momentum-dependent 
effective mass, corresponding to the finite value for the chiral condensate and manifesting the spontaneous breakdown 
of chiral symmetry (SBCS). Moreover, the momentum-dependent quark mass plays the role of a UV regulator by 
construction. Since we are interested in the case at finite T, the caloron with trivial holonomy is used to modify 
the instanton parameters [26] [27]. By doing that, the instanton size p becomes a smoothly decreasing function of 
T, signaling weakening nonperturbative effects of QCD. In addition, we utilize the fermionic Matsubara formula for 
inclusion of T in the relevant physical quantities in hand. In order to take into account the external magnetic field, 
which is produced in the peripheral heavy- ion collisions, the Schwinger method is employed [28H30] . 

Using the Kubo formula for the electrical conductivity, we compute a as a function of T and the relaxation time 
r for the SU(2) light-flavor sector in the chiral limit, i.e. m u ~ rrid ~ 0. From the numerical results, it turns out 
that a « (0.02 - 0.15) fm -1 for T = (0 - 400) MeV with the relaxation time r = (0.3 - 0.9) fm. In addition, 
we also parameterize the electrical conductivity as cr/T « (0.46, 0.77, 1.08, 1.39) Cem for r = (0.3,0.5,0.7,0.9) fm, 
respectively. These results are well compatible with other theoretical estimations, including those from the lattice 
QCD simulations. It also turns out that the external magnetic field plays only a minor role for a even for the very 
strong magnetic field Bq ~ rn^ x 10 and becomes relatively effective for T < 200 MeV. Moreover, we also compute 
the soft photon emission rate from the quark-gluon plasma, using the electrical conductivity calculated. 

The present work is organized as follows: In Section II, we briefly introduce the present theoretical framework. 
The numerical results and related discussions will be given in Section III. The last Section is devoted to summary, 
conclusion, and future perspectives. 



II. THEORETICAL FRAMEWORK 



Now, we make a brief explanation for the theoretical framework for computing the electrical conductivity. We 
note that all the calculations are performed in flavor SU(2) in the chiral limit (m n ,m^) — >> 0. First, the electrical 
conductivity can be defined in Euclidean space from the Kubo formula as follows [14] : 

<V(P) = "E^f / -0^^cAS{k)^S{k + p) lv \ A . (1) 

Here, ef stands for the electrical charge for a light-flavor (/) quark, i.e. (e n ,e^) = (2e/3, — e/3), where e denotes the 
unit charge yj 47tc^em with the fine structure constant «em ~ 1/137. w p indicates the Matsubara frequency for the 
momentum p for a, being proportional to 2ttT. Tr Cj7 is assigned as the trace over the color and Lorentz indices. The 
subscript A at [• • • ]a means that we have introduced the externally induced electromagnetic (EM) vector field A^. 
Note that, in Eq. (IT]), we have taken into account only the connected-diagram contribution, since the disconnected- 
diagram one oc (e u + e^) 2 contributes to the electrical conductivity negligibly [18]. The longitudinal and transverse 
components for the electrical conductivity are defined with Eq. ([!]) as 

a± = lim lim [cTn 22 (p)], ^11= lim lim[<7 33 (p)]. (2) 

In order to evaluate Eq. ([!]) , we make use of the dilute instanton-liquid model [24j [25] . This model is based on the 
idea that the quarks interact with the nontrivial QCD-vacuum configuration, i.e. instanton ensemble, by flipping their 
chiralities. As a result, the quark acquires the dynamically generated effective mass, and spontaneously breakdown 
of chiral symmetry (SBCS) is explained. The model is also characterized by the two phenomenological parameters 
for the light-flavor sectors: average inter- (anti) instanton distance R and average (anti) instanton size p. For vacuum, 
these values are phenomenologically chosen to be (R,p) ~ (1, 1/3) fm [24l [25] . In Euclidean space, the effective chiral 
action (EChA) of the model reads: 

Seff = -S Pc , M In [ip - iM(D 2 )] , (3) 

where Sp c j 7 stands for the functional trace over the color, flavor, and Lorentz indices, while M(D 2 ) for the 
momentum-dependent effective quark mass from the quark zero- mode solution in the instanton vacuum [25 . Note 
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that D indicates the U(l) covariant derivative iD^ = id^ ^e^A^. In this way, we induce the external electromagnetic 
(EM) field to EChA: the Schwinger method [28-30 . From EChA, one can derive the light-quark propagator under 
the external EM field in the momentum space as [31] 



S{K) 



iM 2 K fi + i[M k + \M k (<j-F)\ 



K 2 



Ml 



k 2 + M 2 



(4) 



Here, = + e/A^ and a ■ F = o ilv F iLV '. F^ stands for the photon field strength tensor. The effective quark 
mass and its derivative M k in Eq. Q are denned by [3"T] : 



M k =M 



2 + p 2 k 2 



M k 



8M p 2 



(2 + p 2 fc 2 ) 3 ' 



(5) 



where Mq and p denote constituent-quark mass at zero virtuality and average (anti)instanton size, respectively. From 
Ref. [32 , Mo is determined as about 350 MeV for 1/p « 600 MeV, to reproduce the pion weak-decay constant F n . 

Now, we want to explain briefly how to modify p and Mq as a function of T, using the caloron solution. Details 
can be found in Ref. [27] . An instanton distribution function for arbitrary N c and Nf can be written with a Gaussian 
suppression factor as a function of T and an arbitrary instanton size p for pure-glue QCD [26] : 



d(p, T) = C Nc A^ s ^ p b ~ 5 exp [-(A N T 2 + hnp 2 )p 2 } 



(6) 



We note that the CP-invariant vacuum was taken into account in Eq. (J6|, and we assumed the same analytical form 
of the distribution function for both the instanton and anti-instanton. Note that the instanton packing fraction 
n = 1/R A and p have been taken into account as functions of T implicitly. We also assigned the constant factor in 
the right-hand-side of the above equation as C for simplicity. The abbreviated notations are also given as: 



$ = -61n[A RS p cut ], p = -bln[A RS (R)}, C Nc = 



4.60 e ~ 1 - 68a Rs Nc 



1 
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ir 2 (N c -2)\(N c 
_ 11N C - 2N f 



1)!' 



(7) 



Note that we defined the one- loop inverse charge f3 and /3 at certain phenomenological cutoff p cut and (R) « R. Ars 
denotes a scale, depending on a renormalization scheme, whereas V3 for the three-dimensional volume. Using the 
instanton distribution function in Eq. ([6|, we can compute the average value of the instanton size p 2 straightforwardly 
as follows 33 : 



p 2 (T) 



Jdpp 2 d(p,T) _ [A 2 N T*- 



- 4vf3jn] 



A N T 2 



fdpd(p,T) 2/^n 
where v = (b — 4)/2. It can be easily shown that Eq. ([8| satisfies the following asymptotic behaviors [33] : 



lim p 2 (T) = 

T^0 



^— , lim p 2 (T) = 



(8) 



(9) 



Here, the second relation of Eq. ([9| indicates a correct scale-temperature behavior at high T, i.e., 1/p « A ex T. 
Substituting Eq. ([8| into Eq. (|6), the caloron distribution function can be evaluated further: 



d(p, T)=C p b ~ b exp [-HT)P 2 \ , H?) = ~A N T 2 



A 2 T 4 



The instanton packing fraction n can be computed self-consistently, using the following equation: 

n^(T) = [CI»P, 



(10) 



(11) 



where we replaced NT/V3 — >• n, and T£y) stands for the T-function with an argument v. Note that C and j3 can be 
determined easily using Eqs. (I8J) and (11 ), incorporating the vacuum values for n « (200 MeV) 4 and p w (600 MeV)" 1 : 



FIG. 1: Left: dzv(p, T) on the p-T plane. The trajectory on the curve represents p(T) in Eq. Right: Momentum-temperature 
dependent effective quark mass in Eq. (14). 



C « 9.81 x 10 4 and /3 « 9.19. Finally, in order for estimating the T-dependence of M , one needs to consider the 
normalized distribution function, defined as follows: 



d N (p : T) = 



d(p,T) 
fdpd(p,T) 



Jb—b -rv 



T»(T )exp[-T(T)p 2 ] 



(12) 



Here, the subscript N denotes the normalized distribution. For brevity, we want to employ the large- N c limit to 
simplify the expression for djy(p,T). In this limit, as understood from Eq. (12), dx{p,T) can be approximated as a 
(5-function: 



lim d N (p,T) = 5[p-p(T)}. 

N c ^OO 



(13) 



The trajectory of this <5-function projected on the p-T plane is depicted in the left panel of Figure[T]by the thick-solid 
line. Using Eq. (13), we can write finally the (|fc|, T)-dependent M as follows: 



M k =M 



_2 + k 2 p 2 {T)_ 



(14) 



141). We show the moment um- 
As shown there, the effective 



where we used M ~ 350 MeV at T — ^ 0. Note that p and n are functions of T in Eq. 
temperature dependent effective quark mass in Eq. ( 14 ) in the right panel of Figure [l 
quark mass is a smoothly decreasing function of momentum as well as temperature as expected. 

Considering all the ingredients discussed above and plugging Eq. Q into Eq. 0, we obtain the following expression 
for the transverse components of a with (/i, v) = (1,1) or (2,2) and F12 = Bq in z direction: 



f 



4e 2 f N c 5^ f d 4 k 



I 



(2tt)4 



fc(fc+p) 

2 



Ml 



2Ml(e f B ) 



[k 2 + Ml\[{k + p)2 + M 2 ] 



(15) 



For the longitudinal one (/i, v) = (3,3), one removes the term proportional to Bq in Eq. (15). Since we are interested 
in the finite T case, we introduce the fermionic Matsubara formula for Eq. (15), then have the following expression: 



E 



8e}N c T6^ 



E 



(2^)3 



k +w k (w k +w p ) 
2 



Ml 



2Ml(e f B )< 



[k 2 



wt 



Ml\[k 2 + (w k + w p ) 2 + M 2 } 



(16) 



Here, n > for representing the retarded Green function. For definiteness, we made use of p = (0, w p ) in deriving 
Eq. (16). Note that we have simplified the mass functions in Eq. (16) by replacing the four momentum into three 
one, k — » fc, in comparison to Eq. (15): 



Mfe « M (T) F%, F k = 



2 + p 2 (T) fc 2 ' 



Ml 



8M (T)p 2 (T) 
[2 + p 2 (T) fe 2 ] 3 



(17) 
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where M (T) and p{T) are given in Eqs. (14) and (13). We verified that this simplification does not make considerable 



changes in the numerical results, and make the analytical and numerical calculations much easier. Note that the 
fermionic Matsubara frequency for the loop, is defined with the finite relaxation time r for computing the finite 
electrical conductivity as in Ref. [T4] : 



1 



2t\w„ 



(18) 



Here, w n — (2n + 1)ttT and n e 1. We want to make an assumption: It is clear that the term l/(2\w n \) in Eq. (18) 
becomes its maximum and most effective when n = 0. Hence, it is reasonable to replace \w. 
parentheses in the right-hand-side of Eq. (18). Moreover, if we consider w. 
further as follows: 



| wo | = in the 



1 



1 



2ttTt 



1 



1 

W p T 



27rT, Wk in Eq. (18) can be simplified 

(19) 



In what follows, we will use Eq. (19) rather than Eq. (18) for the Matsubara frequency in the calculations. 



The summations over n in Eq. ( 16 ) can be developed neatly as follows: 

T 1 



E 

n 

E 



[v 2 K + El] 

T 



„ -tanh f — | 
4 V E k \2 V TJ 



[i-M l + ElY (2,/'/ ) 



sech' 



?7Tsinh 



Eh 
rfT 



Eh 



(20) 



where we have defined a notation as 77 = 1 + l/(w p r). Combining Eqs. ([2|, (19), and (20), we have the following final 
expression for a for w p — >• 0: 



E 

/ 

T7T 

T 



e)N c r< 



I 



d 3 k 

J2nf 



tanh {nrEk) 



d 3 k 
(2tt) 3 



sech 2 {tit Eh) 



sinh (27rrEh) 
2ttt 



Eh 



Mi-4Mi(e f B ) 



(21) 



The energy of the quark is given by Eh = (k 2 +M|) 1//2 . We note that we inserted F^(k 2 ) in the integrals in Eq. (21) 
to tame the UV divergence smoothly in integrating over fc, instead of setting a three-dimensional cutoff, cry can be 
easily obtained by putting Bq = in Eq. ( |2T| ). 

Finally in this Section, we want to discuss briefly about the chiral phase transition of the present model. One can 
explore the chiral structure via the chiral (quark) condensate (qq) within the present model [27} [34] : 



(qq) 



-N r N- 



Ml 



(2^) 3 E k 



1 — exp (- 



T , 



1 + exp (- 



T , 



(22) 



From the numerical results given in Ref. [34], the chiral condensate turns out to be a smoothly decreasing function 
of T which is rather different from the universal class pattern of the chiral phase transition, i.e. there appears the 
second-order phase transition for the SU(2) light-flavor sector in the chiral limit m q —> 0. We note that, if the 
meson-loop correction (MLC), which relates to the 1/N C correction, is taken into account, the correct second-oder 
phase transition was observed for the vanishing current-quark mass, giving T c w 160 MeV [34]. Since we have not 
taken into account the MLC contribution in the present work, the chiral phase transition structure of the present 
theoretical framework is distinctive from that of the reality so that the inclusion of MLC can change the numerical 
results for a to a certain extent. We, however, note that the MLC contribution is ineffective for VCC ~ (J^Ju) in the 
chiral limit in general [35 . Hence, it is rather safe to ignore the MLC contribution for the electrical conductivity as 
in the present work, despite of the problematic issue on the chiral structure at finite T as discussed above. We would 
like to leave this interesting task for the inclusion of the MLC contribution for the future. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



In this Section, we present the numerical results for a and make comparisons with other theoretical results mainly 
from the LQCD simulation data. In the left panel of Figure ^1 we show the numerical results of a± (thick) and ay 
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T = 


T 100 MeV 


T 200 MeV 


T 300 MeV 


T 400 MeV 


r = 0.3 fm 


0.020 


0.021 


0.024 


0.031 


0.049 


r 0.5 fm 


0.034 


0.036 


0.040 


0.053 


0.083 


r 0.7 fm 


0.048 


0.050 


0.056 


0.074 


0.116 


r 0.9 fm 


0.062 


0.064 


0.072 


0.095 


0.149 



TABLE I: Typical values of a in Eq. (21 ) [fm -1 ] at B = for different T and r values. 



(thin) for different r values, r = (0.1,0.3,0.5,0.9) fm in (solid, dot, dash, dot-dash) lines, respectively, as functions of 
T. As for the external magnetic filed, we set it to be Bq = x 10 as a trial. Note that this value of Bq is still far 
stronger than that observed in the RHIC experiment [TH3]. From the figure, it turns out that a is a rapidly increasing 
function of T and show the obvious increases beyond T « 200 MeV. By comparing those cases with and without 5 , 
we also find that the effect from the external magnetic field is considerably small and only relatively effective in the 
low-T region T < 200 MeV. Hence, we can conclude from the present model that a± « a\\ = a. For reference, we list 
a values for some typical temperatures in Table [l[ in which one can easily see that a is almost linear for T < 150 MeV, 
and increases monotonically after that temperature. For instance, at T c « 180 MeV, which is close to the transition 
temperature of QCD, calculated from the full LQCD simulations pll37] . we have a = (0.023, 0.039, 0.054, 0.070) fm -1 
for r = (0.1,0.3,0.5,0.9) fm. It is worth mentioning that, in the recent LQCD simulations [38-41 , the transition 
temperature was obtained as T c & 155 MeV, which is lower than those from Refs. [36, 37 . At T c « 155 MeV for 
instance, we have a = (0.022, 0.037, 0.052, 0.067) fm -1 for r = (0.1,0.3,0.5,0.9) fm. Hence, we conclude that there is 
only small change in a (~ afew%) for different T c values below 200 MeV as shown in the left panel of Figure |5J 

For practical applications as in the LQCD simulations [22j[42], it is quite convenient to parameterize the present 
numerical results for a as follows: 

<t(T) = C EU E C mT m , jr~~T £ H- (23) 

rn=l 

Here, we define Cem = XI/ e / ~ 0.051 for the SU(2) light-flavor sector. Obtained coefficients up to m = 3 are given 
in Table [TTJ As shown in the coefficients, the electrical conductivity can be parameterized almost linearly with respect 
to T, i.e. |C 2 , 3 | ~ 0. Hence, we can approximate them as a/T « (0.46,0.77, 1.08, 1.39) C EM for r = (0.3,0.5,0.7,0.9) 
fm, respectively. 

In Refs. [22j and [42 , employing the SU(3) quenched LQCD simulations, it was estimated that a/T = 7Cem for 
1.5 < T/T c < 3 and a/T = (0.4 ± 0.1)C E m for T/T c « 1.5, respectively. Note that there is one order difference 
between these a values, although the temperature ranges are not overlapped. In the left panel of Figure [2] we depict 
these two LQCD values from Ref. [22 (square) and Ref. [42] (circle), using T c w 180 MeV as a trial, although the 
transition temperatures are slightly higher than this value in general in the quenched LQCD simulations. It is shown 
that the data point from Ref. [42] is well consistent with our numerical result for r ~ 0.3 fm. On the contrary, the 
data point from Ref. [22] for T = 270 MeV is much larger than ours for r = (0.3 ~ 0.9) fm. We verified that, in 
order to reproduce it, r becomes about 5 fm in our model calculation as shown in the left panel of Figure [2] in the 
dot-dash line. In Ref. [43J, the characteristic r was estimated using Ref. [22 , resulting in about (2.2 T/T c ) fm with 
a conservative estimate of the QGP medium size. Taking T w 270 MeV = T c x 1.5 MeV, it is given that r ~ 5 
fm, which is in good consistency with our model estimation as shown in the left panel of Figure [2] Similarly, at 
T w 1.45 T c , it was suggested that a/T w (1/3 ~ 1)Cem m Ref- [18 - If we choose Tr~ 180 MeV again, this result 
provides that a = (0.022 ~ 0.067) fm -1 , which is drawn in the left panel of Figure (triangle) and it corresponds 
to r « (0.3 ~ 0.7) fm in comparison with our results. Moreover, the typical time scale of r was given by tT = 0.5, 
giving r ~ 0.38 fm at T « 1.45 T c [18]. Interestingly, this time scale is quite compatible with ours. In those LQCD 
data, it should be noticed that the unrenormalized vector currents for the Kubo formula ~ (J^Jv) were employed in 
Refs. j22[[42], whereas Ref. [18] made use of the renormalized one. Hence, one may expect systematic errors for a 
given in Refs. [22j[42]. To estimate those errors are, however, clearly beyond our scope. 





r 0.3 fm 


r 0.5 fm 


r 0.7 fm 


r 0.9 fm 


d 


0.46 


0.77 


1.08 


1.39 


C 2 [fm] 


4.00 x 10" 6 


6.66 x 10" 6 


9.33 x 10" 6 


1.20 x 10" 6 


C 3 [fm 2 ] 


-4.87 x 10" 5 


-4.87 x 10" 6 


-4.88 x 10" 5 


-4.88 x 10" 5 



TABLE II: The coefficients Ci j2 ,3 for Eq. (23) for different r values. 



7 




FIG. 2: (Color online) Left: electrical conductivity a as functions of temperature T for different relaxation times r = (0.3 ~ 0.9) 
fm. The thick line indicate the case with Bo = (cry), whereas the thin ones for Bo = x 10 (<t_l), in Eq. (21). SU(iV c ) 
LQCD estimations are taken from Refs. [22 (Gupta), 42] (Aarts), [18] (Ding), and [19] (Buividovich) . Right: do ft photon 
emission rate in Eq. (24), 1Z 1 in Eq. (24) as functions of T for different r's in the same manner with the left panel. For all the 
curves in the left and right panels, we have chosen T c « 180 MeV as a trial. 



In Ref. [19], the quenched SU(2) LQCD was performed, and the electrical conductivity was also estimated as 
a = (0.076 ± 0.010) fm -1 at T = 350 MeV with the transition temperature ~ 313 MeV, taking into account T = 
1.12 T c . To be consistent with others using T c = 180 MeV as above, we depict the data point of Ref. [19] at 
T = 1.12 x 180 MeV « 202 MeV in the left panel of Figure [5] (diamond) , although it was evaluated at T = 350 MeV. 
Being different from other LQCD data, Ref. [19] presented the longitudinal and transverse components of a separately 
in the presence of the external magnetic field. Those LQCD data showed that a_\_ « cry beyond T c for arbitrary values 
of £?o, whereas <j± a \\ at T = and the difference between them is enhanced by increasing Bo. Qualitatively, this 
observation of the LQCD results are consistent with ours as indicated by the thick and thin lines in the left panel of 
Figure [2j In our calculations, cry (thin) is smaller than <j^, mainly due to that the negative sign in front of the term 
cx (efBo) 2 in Eq. (|2l|) in the vicinity of T w 0. On the contrary, ay becomes larger than a± at T — in Ref. [T9] . 

However, we notetnat, if we consider the finite number- difference between the instanton and anti-instanton numbers, 
i.e. AN = Nj — Nj ^ [31], which are not not taken into account in the present work and corresponds to Ua(1) 
QCD anomaly via the finite chiral quark density related to CME [SHE], the conductivity increases with respect to the 
eternal magnetic field. In other words, in order to observe sizable values for the electrical conductivity enhanced by 
the external magnetic fields, it is necessary to have Ua(1) QCD anomaly explicitly in the model, in addition to the 
anomaly-independent terms [3T] . 

Beside the LQCD data, there were several theoretical estimations for the electrical conductivity via effective ap- 
proaches using the Green- function method [14] and ChPT [44] . In Ref. [14] , the electrical conductivity was computed 
for finite temperature and quark density, T = 100 MeV and \i = 400 MeV, which may correspond to the future 
heavy-ion collision facilities, such as Facility for Antiprotons and Ions Research (FAIR) of GSI and Nuclotron-based 
Ion Collider fAcility (NICA) of JINR as noted there. By choosing r = 0.9 fm, it was given that a ~ 0.04 fm -1 . Note 
that this value corresponds to r = (0.5 ~ 0.7) in our results for zero density. In other words, by increasing the quark 
density, the electrical conductivity decreases at a certain temperature, as expected. 

Finally, we would like to estimate the (differential) soft photon (w p — >• 0) emission rate from QGP for the dilepton 
decay rates which is related to the electrical conductivity as follows [I8j [22] : 

dR^, 3c^em 
" 7 = ■? u lo™ p &\p\ = 



11, = U m w p ^ = ^aT. (24) 



The numerical results for 1Z 1 in Eq. ( 24 ) are given in the right panel of Figure [2] for different r values as in the 
left panel of Figure [2] We also depict the other theoretical estimations there, using Eq. (24). It turns out that the 
value of 7Z 1 increases rapidly in the vicinity of T = 0. Beyond that, the slope of 1Z 1 becomes relatively stable with 
respect to T. Although one can easily have numerical values for 7Z 1 from the numbers given in Table [T] by multiplying 



3a EM T/(27H), for convenience for readers, 7Z 1 for some temperatures are listed in Table |III| 
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T 100 MeV 


T 200 MeV 


T 300 MeV 


T 400 MeV 


T 


0.3 fm 


1.181 x 10" 


-5 


2.700 x 10" 5 


5.229 x 10" 5 


11.021 x 10" 5 


T 


0.5 fm 


2.024 x 10" 


-5 


4.500 x 10" 5 


8.940 x 10" 5 


18.667 x 10~ 5 


T 


0.7 fm 


2.811 x 10" 


-5 


6.297 x 10" 5 


12.482 x 10" 5 


26.089 x 10" 5 


T 


0.9 fm 


3.600 x 10" 


-5 


8.100 x 10" 5 


16.025 x 10" 5 


33.511 x 10" 5 



TABLE III: Typical values of in Eq. (24) [fm -2 ] at B = for different T and r values. 



IV. SUMMARY AND CONCLUSION 



We have investigated the electrical conductivity of quark matter at finite temperature and zero quark density 
(T ^ 0, fi = 0), using the dilute instanton-liquid model. The instanton parameters, i.e. the average instanton size and 
inter- (instanton) distance, were modified as functions of T according to the caloron distribution. Employing the Kubo 
formula, we computed it in the presence of the external magnetic field in the order of Bq ~ m% x 10. We compared 
our results with the LQCD data as well as effective theories. Important observations of the present work are listed as 
follows: 

• The electrical conductivity is an increasing function of T and the values for a depends on the relaxation time 
r of the quark matter. Note that we modify the effective quark mass into a decreasing function of T as well 
as |fe|, playing the role of a natural UV regulator and signaling the reduction of the nonperturbative QCD 
effects at finite T. Typically, we have a w (0.02 - 0.15) fm -1 for T = (0 - 400) MeV with the relaxation time 
r = (0.3 - 0.9) fm. 

• Recent theoretical estimations including the quenched LQCD data are well reproduced for r = (0.3 ~ 0.7) fm 
for a wide T range. From the numerical calculations, it turns out that the one-order larger <r, estimated in the 
early LQCD data [22 , corresponds to r « 5 fm, which is much larger than other values. This value r « 5 fm is 
well consistent with a conservative estimate of the QGP-medium size. 

• The external magnetic field provides only negligible effects on the electrical conductivity even for the very strong 
one in the order of Bq ~ m% x 10. However, if we consider the Ua(1) anomaly, which is the main source for CME 
and produced from the number difference between the instanton and anti-instanton in the present theoretical 
framework, this situation can be altered to a certain extent. 

• Using the numerical results obtained, the electrical conductivity is parameterized for T = (0 ~ 400) MeV with 
& — Cem(CiT + C2T 2 + C3T 3 + • • • ) and, in this parameterization, the coefficients for T 2 and T 3 are almost 
negligible in comparison to d. As a result, we have a/T « (0.46, 0.77, 1.08, 1.39) C E m for r = (0.3, 0.5, 0.7, 0.9) 
fm, respectively. These results are again well compatible with other theoretical estimations. 

In order to understand QCD at extreme conditions, the transport coefficients for quark matter are very important 
physical quantities to be studied. In the present work, we have shown that the instanton model reproduced compatible 
results for the electrical conductivity with other theoretical calculations including the LQCD data, and estimated it as 
a function of T for different relaxation times of the matter for the first time. Hence, it is quite promising and challenging 
to investigate other coefficients as well within in the same theoretical framework. Other transport coefficients such as 
the shear and bulk viscosities are under investigation, and related works will appear elsewhere. 
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